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Abstract 

There is presently disagreement between theory and experiment as weU 
as between different theoretical calculations concerning the fine-structure 
splitting of the lowest P state of the neutral helium atom. We believe 
that we have found a minor error in the formulas used by Drake et al. 
(Can. J. Phys. 80, 1195 (2002)) in their calculations, and we may have 
an explanation how the error has occurred. To what extent this might 
resolve (part of) the discrepancy is not known at present. 

1 Introduction 

The fine structure of the lowest P state of neutral helium is of great prin- 
cipal interest, since a comparison between theory and experiment might 
yield an accurate and independent determination of the fine-structure 
constant, a. Unfortunately, various theoretical calculations disagree, and 
there is also a significant discrepancy between theory and experiment, 
when using an accurate value of a, determined in other ways, primarily 
from the g-factor of the free electron. Even with a as a free-running pa- 
rameter, it is not possible to match theory and experiment for the two 
fine-structure separations. 

The most accurate experimental results are obtained by Gabrielse 
et al. PP and by Inguscio et al. [2]. The theoretical calculations have 
been performed by Drake and coworkers 53 ■ as well as by Pachucki and 
Sapirstein [3]. 

The calculations of Drake et al. are based upon the works of Sucher [5] 
and of Douglas and KroU . We have now discovered that there is most 
likely a minor error in the formulas of Douglas and KroU and the corre- 
sponding formulas of Zhang [7] on which the works of Drake et al. are 
based. This is estimated to lead to a correction of order Ry, which is 
beyond the accuracy of Douglas and KroU but might be relevant for the 
works of Drake and Zhang. We do not know at present the magnitude of 
the effect and to what extent this might resolve some of the discrepan- 
cies, but this is certainly worth investigating. Below we shall review the 
analyses of Douglas-KroU (DK) and Zhang (Z), which are based on the 
early work of Sucher (S) [5], and point out where we believe the mistake 
has been made. 
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2 The analysis of Sucher, Douglas and 
Kroll 

In his thesis Sucher performed a perturbation expansion of the Bethe- 
Salpeter equation (BSE) '8', which with somewhat different notations can 
be expressed 



'^{x,x') ~ J J J j d'*a;id*a;2d'*a;'id'*a::2 

X Go{x,x';x2,X2) (~i)T.*{x2,x'2;xi,x'i)'^(xi,x'i) (1) 

Go is the zeroth-order two-particle Green's function, dressed with aU kinds 
of single-particle self-energy insertions, and E* is the irreducible or proper 
two-particle self energy. The function Gq is a product of two single-particle 
Green's functions, satisfying the relation 

{i—-hijG{x,xo)=iS'^{x-xo) (2) 

which leads to (S 1.6, DK 2.19) 

(i-^ - hi^ ^i-^ ~ h2^^{x,x') = i jj A^x-iA'^x'^T.* {x,x' ■,xi,x'i)^{xi,x'i) (3) 

where /ii^2 are the Dirac single-particle Hamiltonians. 
With the wave function being of the form 

*(a;, x') = *(r, r, x, x') = e^'^"^ *(r, x, x') (4) 

where T = (t + 1')/2 is the average time and t = t — t' is the relative time, 
the BSE can after a Fourier transform be expressed in operator form (S 
1.12, DK 2.26, Z 1) 

.^*(e)=3*(e) (5) 

Here, 

T = [e/2 + e-hi^ {e/2 - e - fta) (6) 



and 



ff*W = ^S*(e)$ (7) 



where (S 1.31, DK 3.7, Z 14) 

^ = j de*(e) (8) 

is the " equal-time'^ function. 

The interaction g can be separated into a Coulomb part and a "re- 
mainder" 

g = Qa + QA (9) 

leading to 

*(e) = (.^-5A)"'5c*(e) (10) 
From Eq. (0) it follows that (S 1.32), (DK 3.8) 

flc*(e) = ^/.<l> (11) 
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where 7c is the Coulomb interaction. 

Integrating Eq. pip over the relative energy e, leads to (S 1.34, DK 
3.10) 

* = i (12) 



27r' 

The inverse of the operator T is essentially a product of two electron 
propagators, yielding 



(13) 



where A++, A-- are two-particle projection operators for doubly positive 
and negative states, respectively. This leads to the equation (S 1.47), (DK 
3.26) 



fti + /i2+(A++-A__)/c + i y ^DT-^gt.{T - gi^Y^Ic 



(14) 

where 

D = E-hi-h2 (15) 
This is the starting point for the further analysis. 

3 Perturbation expansion of the BS equa- 
tion 

The operator on the left-hand side of Eq. (|14p can be separated into a 
no-pair Coulomb operator 

He ^ hi + h2 + A++IcA++ (16a) 

a Coulomb virtual-pair operator 

Hai =A++/c(l-A++)-A__Jc (16b) 

and a relativity and transverse photon operator 

HA2=iJ^DT-^gA{T-gAy'lc (16c) 

Starting from the no-pair approximation 

He *c = Ec *c (17) 
the Brillouin-Wigner expansion yields the energy contribution 

AE ^ E - E,^ (*c|i^ + vrv vrvrv ■ ■ ■ (is) 

where 

r = rQ(s) = -^^ (19) 

is the resolvent. This leads to the expansion terms (S 2.19-21, DK 3.43, 
Z 28) 

A£;<'' = (*c|//a|*c) (20a) 
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= (*c|i?ArHA|*c) (20b) 

AE^-'^ = (*c|J^Ar//AriyA|*c> (20c) 

etc. 

It can be shown that there is no first-order contribution from Hai, 
and the first-order energy contribution then becomes (DK 3.44) 



(21) 



where 

J = 3a(1-.F-^3a)"' (22) 

Of special interest here is one of the second-order contributions (DK 
3.46, note some misprints) 

Asf = {*c|JYAirj/A2l*c) = {*c|/cA__ri j ^Dr-^jr-^h\^c) 

(23) 

It can easily be shown that A__ YD — A — (DK 3.41). Using the relation 
p6ap and with = E — c — hi — h2, we have Ec — H^ = Dc — A++_DcA+_|_, 
and the no-pair equation (|17|) can be written (DK 3.51) 

(Dc - A++/,) = (24) 
Then the second-order correction Ai?^^' (|23p can be expressed 

ASf) = {*,!(/, -Oi j ^jr-ijjr-i^^^^ (25) 

This can be combined with the first-order correction AE^^^ (|2ip . yielding 

(*el(i?+/c-i:>c)i j ^T^^JT-^IA'Hc) = (*c|(/c+A£)i j ^T-^JT-^I,\<^,) 

(26) 

Here, the AS term differs in sign from (DK 3.54) and (Z 37). 

The reason for the discrepancy between our result here and those of 
Douglas and Kroll and of Zhang seems to be that the latter use the relation 

= S1S2 = {Si + S2) + S^^y^ = D-^ (Si + 52) (27) 

where 51,2 are electron propagators, and the identity (DK 3.50a, note 
misprint) 

to transform the first-order equation (|21l) to 

AB'i) = (*,|Z?-^(1-Ai5/Z))i y"^(Si + S2)J(Si+S2)/c|*c) (29) 
and the second-order correction (|25p to 

Asf ' = {*c|(/c - Dc) D-^ i I ^[Si + &) J(5'i -f 5'2)/c|*c) (30) 
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Then they cancel D^^ in the first equation against Dc in the second. 
Zhang approximates D by Dc in all AS'^' and A_E'^^ expressions, which 
leads to the same cancellation. According to Eq. (|28|l . this leads to an 
error of 2AE/ (DcD), which explains the difi'erence. 

The difference is of the form 2AE x Ai5*^' , which is of order Ry. The 
correction AE represents the difference between the full energy and the 
Dirac-Coulomb energy, which contains the instantaneous Breit interaction 
and therefore of order Ry. This interaction does not contribute to the 
fine-structure splitting, and therefore AiS'^' is order a'^ Ry, making the 
correction of order Ry. This does not affect the work of Douglas and 
KroU, who study corrections up to a* Ry, while it is of relevance for Zhang 
and Drake who go one step further. 

4 Equal-time approximation 

Another possible source of the discrepancy between theory and experiment 
might be the fact that the calculations - following Sucher - are based upon 
the so-called equal-time approximation, where the particles are assumed 
to have the same time. This is in contrast to the covariant Bethe-Salpeter 
equation, where the particles have individual times. It is hard to tell what 
size an effect of this kind might have on the results and to what extent it 
can be visible at the present level of accuracy. Most likely, however, this 
will be the case at some level. 
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